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Abstract 

In enlarging the field content of pure Yang-Mills theory to a cutoff dependent matrix 
valued complex scalar field, we construct a vectorial operator, which is by definition 
invariant with respect to the gauge transformation of the Yang-Mills field and with 
respect to a Stueckelberg type gauge transformation of the scalar field. This invariant 
operator converges to the original Yang-Mills field as the cutoff goes to infinity. With 
the help of cutoff functions, we construct with this invariant a regularized action for 
the pure Yang-Mills theory. In order to be able to define both the gauge and scalar 
fields kinetic terms, other invariant terms are added to the action. Since the scalar 
fields flat measure is invariant under the Stueckelberg type gauge transformation, we 
obtain a regularized gauge-invariant path integral for pure Yang-Mills theory that 
is mathematically well defined. Moreover, the regularized Ward-Takahashi identities 
describing the dynamics of the gauge fields are exactly the same as the formal Ward- 
Takahashi identities of the unregularized theory. 

PACS numbers: ll.15.Tk, ll.lO.Lm, ll.lO.Hi, ll.15.Ha 
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1 Introduction 



To understand the dynamic of pure Yang-Mills (YM) theory [T], in all the range of the energy 
scale, one needs first a gauge-invariant regularization in four-dimensional spacetime. In the 
seventies, Wilson has built a gauge-invariant regularization of YM theory by approximating 
the continuous four-dimensional spacetime as a discrete lattice [2]. In this approach gauge 
invariance is obvious, but the continuous symmetries of the physical spacetime, such as 
Lorentz invariance are clearly lost. The continuous YM action is recovered when the lattice 
spacing goes to zero. In this regime the theory is perturbative, and for each physical quantity, 
the result calculated in lattice gauge theory must match the result given in a perturbative 
continuous regularization, such as dimensional regularization [3]. For instance, as the lattice 
spacing goes to zero, some integrals which arise from the diagrammatic expansion of the 
action of YM lattice perturbation theory, like the tadpole integrals, become not regularized 
in the infrared (IR) and need an intermediate regularization [U [5]. Therefore, since the 
matching between the discrete and continuous regularization is really not obvious, for the 
study of the evolution of physical quantities with the energy scale, ranging from the IR to the 
ultraviolet (UV) domains, a continuous nonperturbative regularization in four- dimensions 
which will preserve both the spacetime symmetries and the gauge invariance of the theory, 
will be very appealing. 

Apart from dimensional regularization, which is in fact purely perturbative, two differ- 
ent classes of continuous regularization were essentially developed as a starting point for a 
nonperturbative approach, since the advent of lattice regularization. l)In the method of 
higher covariant derivative supplemented by the introduction of Pauli-Villars (PV) regulator 
fields [HE], which are necessary to regularize the one loop diagrams, it is the action which is 
regularized. Because of Gribov [7] ambiguities, in this scheme the regularized path integral 
has nevertheless only a perturbative meaning. 2) In the method of stochastic regularization 
[8], by construction it is not the YM action which is regularized with the help of covariant 
derivative, but the second order Schwinger-Dyson (DS) equations, which are obtained after 
integrating out the Gaussian noise field. Since in this scheme the quantization occurs with- 
out Faddeev-Popov ghosts, the regularized DS equations can in principle be used to study 
the nonperturbative aspects of YM theory. However in this scheme the renormalization pro- 
gram is not straightforward. In fact, due to the lack of a regularized action, the regularized 
Ward-Takahashi (WT) identities can not be deduced from the onset, and one must assume 
[9j that they are identical to the unregularized ones. 

In order to be able to apply exact renormalization group techniques |T0] , we require that 
the regularization scheme does not work only on Feynman diagrams, or DS equations, but 
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regularize the YM path integral as a whole, explicitly and not only formally^. One method 
to regularize the path integral of a given theory, is to add a cutoff function in the quadratic 
term of the action, in order that the propagator becomes a rapidly decreasing function of 
the square of the momenta in momentum space. This idea is one of the cornerstones of 
the construction of regularized path integrals which are suited for the derivation of exact 
renormalization group techniques [H]. In this approach, whereas the path integral that 
quantizises the theory is indeed well defined, the regularization breaks explicitly the internal 
local symmetries of the original action. This is because the gauge transformation mixes 
different scales of the gauge field momenta. Only recently was it shown that this drawback 
can be circumvented by enlarging the fields content of the theory : The higher covariant 
derivative regularization is applied to a spontaneously broken supergroup which embeds the 
original physical non-Abelian gauge group [12] . The original YM theory is then regularized 
by higher covariant derivatives supplemented by a finite set of PV regulator fields, if some 
preregularization is used [13]. But, because the PV regulator fields introduced in the action 
have negative norm by definition [13] , the path integral of the theory is only defined formally 
and perturbatively. 

This paper shows that the path integral of SU{3) YM theory can be regularized as a 
whole, in a continuous and gauge-invariant way, by enlarging the field content of the theory 
by a finite set of complex scalar fields, smeared with cutoff functions. In order to do so, 
we will follow the lines of reasoning which have allowed us to regularize the path integral of 
QED [\3\ in a gauge-invariant manner. The action is constructed in its essence with a gauge- 
invariant building block Ajj^^. This vectorial operator of dimension one is a function of the 
gauge field and of the smeared complex scalar fields, and converges to the standard gauge 
field A^^ when the cutoff goes to infinity. The variation of this vectorial operator under a 
local gauge transformation of the gauge field is compensated by a group transformation of the 
smeared complex scalar field which depends on the cutoff scale and which plays a role similar 
to that of Stueckelberg compensating fields [16]. The regularization is implemented in the 
action by smearing these invariants by smooth cutoff functions. Then all interaction vertices 
of the regularized action are smeared with a cutoff function and the quadratic terms in the 
fields can be inverted to give rise to gauge and scalar fields propagators which are rapidly 
decreasing functions of the square of the momenta. Since the measure and the regularized 
action is invariant under both the gauge transformations of the vectorial and scalar fields, 
when the integration over the scalar fields is being assumed implicitly, the regularized WT 
identities relative to the gauge field are exactly the same as those which can be deduced 
naively from the unregularized action without gauge fixing. 

The paper is organized as follows. Section II is devoted to the explicit construction of the 

^For instance in the case of dimensional regularization the path integral has only a perturbative meaning. 
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regularized gauge-invariant action of the YM theory with the help of the invariant vectorial 
operator Aj^"". The flat measure relative to the scalar fields is showed to be invariant under 
the gauge Stueckelberg type transformation. In Section III, in order to be able to invert 
the gauge vector kinetic term without breaking explicitly the whole gauge invariance of 
the regularized action, we add to the action a gauge-invariant part which is built with the 
invariant operator Ajj^^ and which is reminiscent of the standard covariant gauge fixing 
term. Here the addition of a mass term for the real and imaginary parts of the complex 
scalar fields is required, so that the full action remains regularized. In Section IV, we show 
that the mass term relative to the real part of the complex scalar field, that is not invariant 
under the Stueckelberg type of gauge transformation, is in fact harmless. This implies that, 
after integration on the scalar fields, the regularized WT identities describing the dynamics 
of the YM gauge fields are exactly the same as the formal WT identities of the unregularized 
theory without gauge fixing. 

2 The regularized gauge- invariant action 

If = A'^T", is the SU{?>) color gauge vector boson, we define the 
invariant operator 

with the generators T° related to the Gell-Mann matrices by 

= i^. (2.2) 

Here, ?7 is a 3 x 3 dimensionless complex matrix acting on the generators T", and hereafter 
an expression like d^AB means that the partial derivative acts only on the factor A. By 
construction the vector A]!^"" is anti-Hermitic and invariant under the local transformation 
of A^ and U by an arbitrary SU (3) matrix 

A^ VA^V+ + d^VV+ (2.3) 

U VU. (2.4) 

The invariance of A^^^ relies only on the property of unitarity of the matrix V . Thus, the 
unitarity ^ of the matrix U is in general not required. In order to construct a regularized 
gauge-invariant action, we start with the gauge-invariant expression 

^;rW = j dx {p^^'\z,x) {d.A'r - d.A';^) -p2{z,x)[A';\A'r]]. (2.5) 

^liU were unitary, the definition (|2.ip is that of the gauge- invariant operator given in T7]. 



matrix valued gauge- 
(2.1) 
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Here and hereafter we work in Minkowski space. We choose the signature of the metric to 
be (1, —1, —1, —1) and the notation dx = d'^x and dk = d'^k/{2'n'Y. In the expression (12.51) 
the ultraviolet UV or infrared IR regularization are implemented through the real scalar 
functions 

Pi{x,y) = [ dke-''^'-y^p,{^,K) (2.6) 



pr'ix,y) = / dk , (2.7) 

K and A being respectively the IR and UV cutoff scales. The term proportional to the 
function p[ will contribute to the kinetic part of the action, and hence to the inverse of 
the propagator. As a result, when the path integral of the theory is expressed in momentum 
space, the infinitesimal real parameter t], which was only introduced to give a mathematical 
meaning to the Fourier transform (12. 7p . can be set to zero . In Euclidean momentum space, 
the cutoff functions Pi{k) are rapidly decreasing functions of k"^ in both the IR and UV 
domains and verify the condition ^ 

lim pi{k,A,K) = 1. (2.8) 

A^oo,K^O 

By definition, the functions (12. 6p are regularized forms of Dirac's S function, i.e. 

lim p^ix,y) = 6{x-y). (2.9) 

The expression (12.51) is apparently quartic in the matrix U. In order to simplify and to 
express the operator (12.51) in terms of the usual YM stress tensor, we suppose as a first step, 
that the matrix U is unitary. Then using the relation 

d^U+U = -U+d^U (2.10) 

(12. 5p becomes 

F^riz) = JdxS^U+F;^S'i^^x)U+ip[''\z,x)-p2{z,x))[{d,U+AM + U+A,d,U 

-d,,U+dM)-ip^u)]\, (2.11) 



where the regularized YM stress tensor is defined by the expression 

(z, x) = p[-'^ {z, x) {d^A, - d,A^) - p2{z, x) [A^, A,]. (2.12) 



•^Notice that there are a large number of cutoff functions which arc valid in the IR and UV domains 
and which converge to the Dirac's S function as the cutoff A goes to infinity, whatever the value of k. For 

,2 Jl 

instance, this is the case for the cutoff function — - — . 
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If one relaxes the condition of unitarity of the matrix U, one can check that the new ex- 
pression (12.111) of the operator (12. 5p remains invariant under both transformations (12. 3p and 
(12. 4p . Notice that in this case the relation (I2.10p does not hold true, and hence the rela- 
tions (12. 5p and (12. lip define two distinct invariants. We select the last expression (12. lip 
as the building block of the regularized YM action. In the limit A +oo , k — » 0, one 
can set the parameter t] to zero, and the right-hand side of (12. lip converges to the operator 
U+F^M- The trace of the operator f/^F^"''^'^ does converge to the standard, unregularized 
Lagrangian density of QCD if the matrix U is unitary. 

The next step consists in choosing the form of the dimensionless matrix U. Since the ma- 
trix U is needed to restore the gauge invariance in (12. lip , it can be expressed as a functional 
of some scalar fields which will play the role of Stueckelberg compensating fields [161 [IT]- 

First to avoid any nonpolynomial interaction between the Stueckelberg fields, which will 
enter the action we are looking for, second to have the path integral measure invariant under 
the transformation of the Stueckelberg fields that induces the transformation (12. 4p . and third, 
so that the invariant operator Ajj^^ converges to the standard gauge vector as the cutoff 
A goes to infinity, we define the 3x3 dimensionless matrix U as the linear combination 

U = I+{^. (2.13) 

Here, the matrix is expressed in terms of the complex scalar fields (pa as 

= 0aT^ (2.14) 

and by convention, we use Latin letters when the index runs from to 8. In the definition 
(I2.14p . except the matrix which is given in term of the identity matrix / as 

r° = (2.15) 

the eight remaining matrices are defined in (12. 2p . With the definition (I2.15p . the matrices 
T"" are normalized as 

TiiT^T") = (2.16) 

and the properties of the matrices (12. 2p are recalled in Appendix A. We do not impose any 
constraint on the nine complex fields 0j, so that the matrix U is not unitary in general^. 
Due to its definition, this matrix is an element of the Lie algebra of the U{3) unitary group. 
Since the matrix U transforms like (12. 4p . the transformed matrix 0' verifies the equation 

(f)' = V(f) + A{V -I). (2.17) 



^For a unitary matrix [/, the invariant measure of the fields 0^ is hard to define, as it is the case in the 
nonhnear sigma model |14j . 
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It is clear that the set of transformations (12.171) , acting on the matrix fl2.14p , is a nonhnear 
representation of the SU{3) group. In fact, if g is the group element of SU{3) represented 
by the matrix V, the transformed field 0' can be written as 

0' = 30, (2.18) 

and we have 

9\9^) = 9'9^^ (2.19) 

with the identity group element e represented by the identity matrix J, and the inverse group 
element represented by the matrix V^^. 

How does the measure of the scalar field transform under the transformation (I2.17P ? 
As show in Appendix B, if $ is the nine- dimensional complex vector of components 0j, the 
transformed fields 0^ are given according to $ by the vectorial relation 

$' = S-^T {V) 5$ + AS-^T {V - I) Ci. (2.20) 

Here the nine- dimensional matrix representation T{V) of the SU{3) gauge group and the 
constant vector Cj are respectively defined in (IB. II) and (1B.2|) . and S is the constant matrix 
defined by the relation (IB. 61) . 

Then the Jacobian of the transformation (12.201) is given by 

= S{x - x') det [S-^T {V) S) = 6{x - x'), (2.21) 

because the determinant of T {V) is unity (IB. 41) . In the same manner, one can see that 
the Jacobian of the transformation of the Hermitic conjugate scalar fields (pf, induced by 
the transformation (12.171) . is also unity. As a result, since the scalar fields 0j and (pj' are 
independent variables the measure 

i=S i=S 

v<f)V<p+ = n ^'^^ n ^'^^ (2.22) 

is invariant under the transformation (j27T7j). 

Since the form of the matrix U is now given by the expression ( I2.13p . the invariant tensor 
( 12. lip is recast in terms of the regularized YM stress tensor (12.120 in the expression 



F^n^) = I dx , 



z, x) + p2 (z, x) Gf,^ + (^p[ {z, x) - p2 {z, x) j if^^l , (2.23) 



where the antisymmetric tensors G^v and H^i, are given in terms of the gluons fields A^, the 
standard YM strength tensor F^y and the complex scalar fields by 

G^, = i(0+F^. + F^.0) + ^0+F,..0 (2.24) 
H^, = [i9^(0+A, + A,0) + Jj9^(0+A,0) - ^9^0+9,0] - [/i ^ z/] . (2.25) 
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A candidate for a regularized form of the YM action can be defined as 

S{A, 0+, 0) = ^ Tr 1 dzF'^^ ^'^{z)Fl:\z), (2.26) 

g being the dimensionless strong couphng constant. The exphcit form of this action can be 
expressed as 

5(A,0+,0) = 5yM + 5A<^+0, (2.27) 

where the pure gauge part of this action is given in terms of the regularized YM stress tensor 
f l2.12p by the expression 

Sym = Tr y dxdydz F^^^ ^^{z, x)F«^^(z, y), (2.28) 

and where Sa4,+4, is given in terms of the operators G^^ f l2.24p . H^^ (12.251) and Fl^^^ (12.121) 
by, 



Sa<i>+<p = dxdydz ^p2{z, x)p2iz,y)G'"'G^i,{y) 

+ {p[''\z,x) - p2{z,x)){p[-'\z,y) - p2{z,y))H'^''H,,{y) 



+2 



p,{z,x)G^''F;^:^{z,y) + {p[''\z,x) - p2{z,x))H'^''F;i^^{z,y) 



+P2{z,x){p[~'\z,y) - p2{z,y))G^''HM] I. (2.29) 



Under what conditions is the action (12.271) regularized ? In order to study this point, we 
rewrite the action (I2.29P as a sum of five terms 

SA<f>+<f> = Sgg + Shh + Sgf + She + Sgh, (2.30) 

each term being respectively associated to the product of the operators which enter the 
expression of the action (12.291) from the left to the right. The Fourier transform of these 
terms are given in Appendix C. First of all, we suppose that the cutoff functions pi and p2 
(12. 6p are chosen in such a way that the product p^Vs are still rapidly decreasing functions 
of k"^ in momentum Euclidean space. The simplest choice is 

P2{k)^pl{k). (2.31) 

Then, since the kinetic term arising from the term F^^^^ (12.280 and (12.120 generates a propaga- 
tor ^ proportional to pi in momentum space, all the terms of the action (12.270 are regularized 
with respect to the gauge field, except the terms of (12.291) which do not contain the factor 



^In the sequel, we suppose that some gauge fixing mechanism will be implemented in order that the kernel 
of the kinetic term can be inverted. 
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P2, i-e. the terms contained in the expressions Shh and Shf (12.301) . The divergent function 
p^^ enters in both the Fourier transform of these expressions, which are given respectively 
in Appendix C by (IC.lll and (1C.2I) . In order to regularize them, we assume that the field 
is in fact a smeared field defined by 

= y" dz p^{x,z)ip{z). (2.32) 

In addition, one requires, 1) that the function ps, which enters in the definition (I2.32p of the 
smeared field 0, behaves like 

p,{k)^p\{k), (2.33) 

as the cutoff A goes to infinity, and 2) that the divergent cutoff function Pi^{k) occurring 
in the expressions ( IC.ip and ( ]C.2[) . must be multiplied by the larger number of the smeared 
field which is a function of the variable k. This means that, in order to calculate the 
Fourier transform of the expressions Shh and Shfi the Dirac delta function that reflects 
the momentum conservation at each vertex, must be integrated over the momentum variable 
of the field first. In that explained in Appendix C, the shift of the momentum 

variable of the field is in general not allowed in some terms of the expressions (1C.1I) and 
(1C.2I) . if these terms become not regularized. 



Keeping in mind the rules outlined before, we see that both the expressions fIC.ip and 



flC.2l) are indeed regularized in the gauge field sector, if the propagator of the gluon field is 
a rapidly decreasing function in Euclidean momentum space. Moreover, given in Appendix 
C, the explicit form in momentum space of all the interaction terms between the gauge field 

and the smeared bosons fields (12.321) . which contribute to the action (12.271) and (I2.30p . 
it is clear that the action (12.271) is now completely regularized in the gauge and scalar fields 
sector as long as the propagators of the ip fields (12.32^ are also rapidly decreasing functions 
in Euclidean momentum space. 

Up to now the action (12.271) does not contain any sector that is purely quadratic in the 
scalar field y?. In order to define the gluon and scalar field propagators, we will add to the 
action a new term, which is invariant under both the transformations (12. 3p and (12.171) of the 
gauge and scalar fields. This invariant term enables us to invert the quadratic part of the 
gauge field A^. Since this new term contains a quadratic part in the field which does not 
enforce the regularization, we will also add to the action (I2.27P a mass term for the fields (/? 
which will render the theory completely regularized. 



3 The scalar field sector 

Our aim is to add to the action (I2.27P a new term in order that the gluon kinetic term can 
be inverted without explicitly breaking the invariance under both the gauge transformations 
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(12. 3p and (I2.17P . The new action, thus constructed, must be gauge-invariant in all the range 
of the cutoff A, regularized and must converge to the standard one of QCD as the cutoff goes 
to infinity. Because the part (12.271) of this new action is already converging to the standard 
form of QCD as the cutoff goes to infinity, the searched new term must vanish as the cutoff 
goes to infinity. Then, the part of the action which will give rise to the covariant form of the 
gauge field propagator can be defined as 



''Gauge 



^Tr / f/xrfyrfz(p^^)(^,x)-p2(^,x))(pS-')(^,y)-p2(^,Z/))5M;™9^M^™(y), (3.1 



the gauge-invariant operator A^™ being defined in (12. ip . the scalar functions Pi~^'* and p2 
being respectively defined in (12.70 and (12. 6p . and ^ is a real parameter. If we recall that the 
operator U which enters in the definition of A^^'" is expressed in terms of the smeared field 
(f) (I2.32P by (I2.13p . we obtain in momentum space ^ 

SGauge = ^ Tt j dk {p^\k) - p2{k)f \k^k,A^'{k)A'{~k) + ^(0+(A;) - <P{-k)) 

x(0+(-fc) -0(A;))1 +5Ga«,e(A0), (3.2) 



where the term SGauge{A.,(t)) is explicitly given by the expression flC.3p . Adding both the 



quadratic part in the gluons fields of the expression (13.20 . to that of the regularized YM 
action Sym (12.280 . allows us to invert the gauge field kinetic term, and thus to define the 
free gluon propagator. This propagator is given in momentum space by 

{Afj_aA^fs) = —pl{k)-^5af3[{k'^g^„ — kfj^ky) + ^(i_p^(fc)p2(fc))2 ^t^^A • i'^-'^) 

What happens in the pure scalar fields sector ? We define the smeared field in terms 
of the smeared fields a and vr by 

= a + i7r, (3.4) 

and, in the same way that the smeared field (12.140 is associated to the field in (I2.32p . 
we define the anti-Hermitic matrices q and w as 

= <j + iw, (3.5) 

whose components q and on the basis formed by Tq and the eight generators Tq,, are real 
scalar fields. Because the fields <,a and Wa are real, we have, 

^+(fc) = -<^(-A;) (3.6) 
ro+(fc) = -vj{-k). (3.7) 

fi • • ■ ■ ■ ( — 1) 

Remember that the hmit 77 ^ must be taken in the Fourier transformed of the function p\ 
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In that case, the second term of (13. 2p will give only a contribution to the propagator of the 
<^ field. Up to now, since there is no other quadratic term in the pure scalar fields sectors 
of the parts of the action given in (12.271) and (13. 2p . and knowing that the cutoff function 
which enters in the definition of the smeared a field (13. 4p behaves like (12.330 . we find that 
the free propagator of the ^ field (13. 5p will behave respectively in the IR and UV domains as 

(^^) - ^Pi\k). (3.8) 

Such a behavior for the q field will then spoil the UV regularization of the theory and worse 
will enforce the IR divergences. Moreover a perturbative calculation cannot be performed 
because in the actions (12.271) and (13. 2p there is no quadratic term in the w fields (13. 5p . 
To cure this disease, we will add to the action a mass term for both q and w fields, while 
preserving if possible the invariance of the action under both the transformations (12.30 and 

dZIZD- 

We will first construct a mass invariant term for the field w which is proportional to the 
sum (f)^ + (f). Since the operator 

f/+f/-/ = ^(0+ + 0) + ^rt, (3.9) 

where U is defined by (I2.13p . is invariant under the transformation (I2.17p . the searched mass 
term will arise from the following action 

5. = -M^Tr j dxdy p^-'\x,y){ct>^ + ct>+lct>+ct>){<P^{y) + <P{y) + {cp+{y)ct>{y)), (3.10) 

where the function p\ is defined in (12.70 and M is a mass parameter. Therefore if we use 
the definitions (12.321) and (13.51) . and the relations (13. 6p and (13.71) . and suppose for simplicity 
that the cutoff function ps (12.330 . which enters in the definition of the smeared fields a and 
TT (13.40 . is in fact given by 

p,{k)=pl{k), (3.11) 
the action (I3.1O0 becomes in momentum space 

= M1tt:S^J dk pf{k)w{k)w{-k) -{j dkdp ^^{pMp + k)w{-k) 

j dkdpdqpl{k)ip+{p)ip{p-k)ip^{q)ip{q + k)^. (3.12) 

In order that the second term of (I3.12p gives rise only to regularized contributions, the 
free propagators of the <^ and w fields must be rapidly decreasing functions in Euclidean 
momentum space. 
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Regarding the <^ field, there exists a priori no mass term invariant under the transforma- 
tion (12.171) . For simphcity, we choose the following term, 

^. = |gTr / dxdy p[~'\x,y){<P+ - <P){<l>+{y) - <P{y)), (3.13) 



4g^ 

for the <^ field mass term, which is given in momentum space by 

Af2 rr„ / JU -2, 



S„ = ^Ti j dk p^\k)^{k)<;{-k). (3.14) 

We will see in the following, that the breaking induced by the term (I3.14p does not really 
spoils the gauge invariance of the theory. Now from the relations (13.21) and (I3.14p . we obtain 
for the free propagator of the <^ field 

{^^) = ^ — • (3.15) 

From the quadratic term of the action (I3.12p . and from the relation (I3.15p . we conclude that 
the free propagators of the w and <^ fields have the requested behavior in the IR and UV 
domain, i.e. 

~ ~ (3.16) 

4 The regularized gauge- invariant path integral 

From the preceding sections, we conclude that the form of the regularized action of the YM 
fields is given in adding to the action (I2.27P both the contributions (13. ip of the gauge sector 
and that of the mass terms (I3.1QP and (13.130 of the w and <^ fields. Then the regularized 
action is 

SR^g{A, (^+, y?. A) = SinM, (^+, LP, A) + (4.1) 

where 

Slnv^-^i V^^, V^, A) = Sym + 5'^(/>+0 + Scauge + S-jt-, (4.2) 

is invariant under the gauge transformations (12.30 and (12.40 . and we recall that A is the cutoff 
scale and M is the mass of the w and <^ fields. In the configuration space the regularized 
path integral which is normalized to unity reads 

ZReg{l A) = /^-^^^^^^^-^ f \ (4.3) 

where S{J) is given as usual in terms of the external sources for the gluons fields by 

S{J) = -ITt I dxJ^A^\ (4.4) 
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By construction, when the cutoff A goes to infinity, the action (14.11) has the property 

^im SR,g{A, ip+, <f,A) = ^TTj dxF^^'^F^, - ^ Tr ^ dx^^^. (4.5) 

This tells us that the auxiliary bosons fields if will decouple from the theory. We will show 
that this decoupling occurs when the mass M of the bosons fields (p are at least greater or 
equal to the cutoff scale A. 

After a loopwise expansion of the path integral (14. 3 p with respect to the part of the 
action which contains the vertices of the fields and integration over these fields, as show in 
Appendix D, the proper vertex V ( ID. 71) associated to any gluons effective amplitude behaves 
like 

T ^ {^y^"'0{K-'), (4.6) 

as all the internal momenta scale with the cutoff A. Here n and I^p are respectively the 
number of pure (p fields vertices and the number of internal (p fields occurring in F. As 
show in Appendix D, the number k flD.81) is a positive integer for all the gluons effective 
vertices, except for the two, three and four legs effective vertices, which are proportional to 
the vertices of the pure Yang- Mills action flC.4p . where k is zero. As a result if 



K A^ < or < A^ < (4.7) 

the auxiliary fields (p decouple from the gauge sector. Notice that the rescaling 

^ A0, (4.8) 

of the smeared field cj) (I2.32p does not change this conclusion. This is because, in each internal 
(p boson line the rescaling of the field is exactly compensated by the rescaling 

M KM, (4.9) 

of the mass M entering in the propagator, so that the behavior of the proper effective vertex 
r (14.60 is left unchanged. In order to have only one scale in the regularized path integral 
(14. Sp . a good choice is to take for the mass M of the ip fields the cutoff scale A. 

If we denote pi by p, and remember that the cutoff function p2 is constrained by the 
relation (I2.3ip . and that we have fixed the function p3 by the equation (13. lip , the action 
(14. ip is still regularized if we choose the following relations between p2 and ps 

P2{k)=p,{k)=p{kf. (4.10) 

What are the WT identities that can be deduced from the expression of the path integral 
(14. 3 p ? The WT identities follow from the invariance of the path integral (14. 3 p under the 
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shift of variable induces by the infinitesimal form of the gauge transformation (12.31) and the 
smeared field transformation (I2.17p . i.e. 

^ - C^^^A^e-, + 9^e« (4.11) 
^ + e0 + Ae, (4.12) 

where we use the shorthand notation e = e^T". By construction, all the terms of the action 
(14. ip . except the term S„, are invariant under both the gauge transformation (12.30 and the 
smeared field transformation (12.171) . In addition the measure Vip^Vip is also invariant under 
the transformations of the fields (p which are induced by the transformations ( 12.170 of the 
smeared field (f) (12.321) . This is because, 1) the measure T>ip~^T>ip is proportional^ to the 
measure and 2) the transformation (12.171) of the smeared field (p (12.321) leaves the 

measure (I2.22p invariant. 

In order to derive the WT identities for the gluons fields, we will show that the symmetry 
breaking term (13.131) is harmless. In fact, if we introduce the auxiliary matrix valued field 

f? = f?an (4.13) 

whose components Qa are real scalar fields, and define the action 

S^, = -^Tt J dy[2giy) J dxp^-'\x,y)a + g'iy)], (4.14) 
the path integral (14. 3 p can be written as 

^i(s^„„(A,<p+,V=,A) + Sa£, + S(J)) 

This is because the Gaussian integration over the g fields transforms the action S^g (14.141) 
in the action (I3.14p . Since under the infinitesimal transformation (I4.12p the components 
of the smeared fields a and vr are mixed in the following way 



rr,A -r, +-n -n i{Si„^(A,v^ ,,p,A)+Sae + 3(J) ] 



S(To = -^TCaea (4.16) 
St^o = (4.17) 

the infinitesimal variation of the action S^^g (I4.14p under the transformation (I4.12p is given 
by 

SS^g = j dxdyp^~^\x,y) -^0oiy)earCa + 0(3{y){leai(^jCaf3'y + T^jd^py) 

+ (^7ro-A)64]. (4.18) 



^This fact is readily seen in momentum space. 
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The invariance under the infinitesimal transformations (14.111) . (I4.16P and (I4.17P of the path 
integral (14.151) gives the following identity 



- liC^P,J'pA,, + d.JH)) = 0. (4.19) 

In a intermediate step, and in order to express the connected Green's functions, we introduce 
respectively in the path integral (14.151) the external sources Ka, La and Na for the scalar 
fields ^a, (13.51) and Qa (14.131) . in replacing the source term S{J) (14.41) by the expression 

5( J, K, L, N) = -f Tr j dx {J^A^" + Kq + Lw + Ng). (4.20) 

Then, the generating functional of connected Green's functions in the presence of the external 
sources for the gauge and scalar fields is defined as 

ZR,g{J, K, L, A^, A) = e^«=«(-^'^'^'^'^). (4.21) 

As usual, we define for each field its vacuum expectation value in the presence of the external 
source. These C-number functions are obtained by taking the functional derivative of W^eg 
with respect to their own sources. We also recall that the generating functional of IPI 
functions is defined in terms of the vacuum expectation values of the fields by the Legendre 
transformation 

TReg{A^ w, g, A) = -iWReM L, iV, A) + ^ Tr ^ dx (J^A^ + K<, + Lw + Ng). (4.22) 

In that case there is a one-to-one correspondence between the classical field thus obtained 
and its external source. For instance these conjugate relations hold true for the fields A^ 
and g 

=-^^^ ^M"= - (4-23) 

We are now able to show that the functional derivative of (-jf^) (14.181) with respect to 
the vacuum expectation value of the gauge field A^ vanishes identically. 

At first, when taking the functional derivative of {-^) (I4.18P with respect to the classical 
field Af^, the terms which are only proportional to (g) do not contribute. This is because, 
1) the vacuum expectation values of the fields A^ and g are by construction independent 
variables, and 2) at the end of the calculation we will set to zero the external sources relative 
to the scalar fields. Then, using the chain rule which expresses the functional derivative ^ 
in term of jj, it remains to show that the functional derivative of the connected Green's 
functions {gaT^b) and {ga^b) with respect to the gauge fields sources vanish effectively. 
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The equation of motion of the auxihary field ^, which is deduced from the expressions of 
the generating functional (14.211) and of the action S^g (I4.14p , can be expressed with the help 
of the relations (I4.24p in terms of the classical fields as 

M'[J dy pix,y)<^M + Qa]-g''-^ = 0. (4.25) 

If we take the functional derivative of both sides of the equation (14.251) . first with respect to 
the classical fields ruhiy) or <jfe(?/), and then with respect to A^{z), we get 



= 



6A'i{z)5mt,{y)Se, 



Because vr and a are respectively the smeared field associated to the scalar field w and ^, the 
vanishing of the IPI functions (I4.26P shows that the functional derivatives of the connected 
Green's functions {QaT^h) and {Qa'^b) with respect to the source J of the gauge field also 
vanish identically, whatever the values of the scalar fields sources. As a result the functional 
derivative of the vacuum expectation value {^jf^) (I4.18P with respect to the classical field 
indeed vanishes. Then, in expressing the WT identity (14.190 in terms of IPI functions, 
we get the following relation 

sl/ivTsAf^^'^isA'^ - C^f3^^^^6{x -y)+ d'^j^^^^ = 0. (4.27) 

In this identity T^eg is the generating functional of IPI functions defined by 

TR,g{A, A) = -tWnegU + j dx J^A^, (4.28) 

which means that the integration over the scalar fields ip has already been done, at least 
implicitly, in the path integral (14. 3p . The regularized WT ( I4.27P is our main result, and shows 
that all regularized IPI functions which are deduced from the regularized path integral (14.31) 
are all transverse with respect to the gauge field when the external source is switched off. 



5 Conclusion and outlook 

We have shown that the path integral of pure YM theory can be regularized in four- 
dimensional physical space in a gauge-invariant manner. The regularization is implemented 
nonperturbatively at the level of the action through cutoff functions. The price to pay, in 
order to maintain the original gauge invariance, is to enlarge the field content of the theory 
by a set of auxiliary complex scalar fields that play the role of Stueckelberg compensating 
fields. By contrast with any PV inspired regularization, the Gaussian metric induced by 
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these fields has the good sign, and the path integral is here a mathematically well-defined 
object. Since the action, and therefore the path integral of the theory, is regularized nonper- 
turbatively in all the energy range, this regularization will be useful to study the evolution 
of physical quantities with the energy scale, without making any assumptions for matching 
of the quantities calculated in the IR and in the UV domains, as is the case in lattice gauge 
theory. In a forthcoming paper, we shall study first the consequences of this regularization 
in perturbation theory, and then beyond. 
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Appendix A 

Here for convenience and to fix our conventions we recall the commutation relations of the 



generators (12. 2p of the group SU (3) 

[T",T^] = C„^^T^ (A.l) 

{T",T^} = -l6^f, + td^p^T\ (A.2) 

the structures constants C^/j-y having the opposite sign of those defined in [18]. We recall 
that the trace properties of these matrices follow recursively from the relation 

T-T" = + liC^p^ + id^p,)T\ (A.3) 

For instance the trace of the product of three and four generators are given by 

Xr(T-T^T^) = -\{C^^, + id^p^) (A.4) 

Tr(T-T'^T^T^) = ^S^^pd^s - + id^p^,){C^sa' + td^So.')- (A.5) 



Appendix B 

To have an explicit expression of the transformation of the scalar fields (pi induced by the 
matrix transformation (12.171) . we recast this equality between two 3x3 complex matrices 
as an equality between two complex vectors having nine components. Let Cy be the nine 
components vector column obtained from any 3x3 matrix V by putting the successive 
columns of the matrix as a single column in increasing order from top to bottom. Then, 
the tensorial product involving the complex 3x3 matrix V 

T{V)=I®V, (B.l) 
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is a 9 X 9 matrix T{V) that relates the constant column vector Ci associated with the 3x3 
identity matrix / to the column vector Cy as, 

Cv = T{V)Ci. (B.2) 

By construction, the linear mapping T between the set of 3 x 3 complex matrices V and the 
set of 9 X 9 complex matrices T{V), is an group isomorphism, that is to say, for any matrices 
V and V the relations 

T{VV') =T{V)T{V'), T{I)=I and T{V+)=T+{V) (B.3) 

hold. Moreover, due to the definition fIB.ip we have the obvious relation 



detT{V) = {detvy. (B.4) 



Using the definition flB.2|) and the properties flB.3|) . the vectorial form of the equation (12.171) 
is then given by the equation 

T{(f)')Ci = T{V)T{(f))Cj + AT{V - I)Ci. (B.5) 

We can write the vector T{(j))Ci associated to the matrix fl2.14p as 

T{(j))Ci = S$, (B.6) 

where S is an invertible 9x9 constant matrix acting on the vector column $, whose com- 
ponents are the fields (pi defined in (12.141) . the index i running from to 8. In that case, we 
can rewrite the equation ( ]B.5|) as 

= T{V)S^ + AT{V - I)Ci. (B.7) 



Appendix C 

In this appendix we give explicitly the terms of the action (I2.27P in momentum space. In all 
the terms, the integration over the Dirac's delta functions which result from the translational 
invariance of the action are first performed with respect to the momenta of the (p field, and 
owing to the conditions (I2.3ip and (12.331) the product of cutoff functions V2 and Pi^ps are 
rapidly decreasing functions of the square of the momenta in Euclidean space. Moreover, as 
said before, the inverse Fourier transform of the scalar function (12. 7p p\~^^ is defined when 
the infinitesimal parameter 77 tends to zero. After a lengthy but straightforward calculation. 
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we obtain respectively for the terms Shh-, Srf and Scauge which are defined in (12.291) , (12.301) 
and flO) . 



S 



HH 



' Tr j dkdpdq {p-,\k) - P2{k)f\{k''q^ - kqk,)A'^{p) - A;)T" 

+r"0(A; - p))0+(g)0(g ~ k) + § j dp'<P+{p')T''<P{p' + k- p)(p+{q)(p{q - A;) 
Hk'g,. - k.kMMA'^iq) - k)T- + T-<P{k - pm+{q + k)T^ 

+T^{-q -k)) + ^ [ dp'{(P+{p - k)T'' + T'^^ik - p))(t)^{p')T^{p' - k - q) 



+ ^ / dp'dq'<P+{p')T''<P{p' + k- p)<p^{q')T^<p{q' - k - q) 



-j4P^{k'^q^ - kqk^)<p^{p)(j){k+p)(i)+{q)(i){q - k) 



(C.l) 



Shf = ^Tr j dpdqp^\p){p^\p)-p,{p))}^^{p\-pqp,)A^^{p)T»<P+{q)<P{q 
Hp'g,. - P,Pu)A>^ip)A''p{q) [i(r"TX-p -q)+ T^T>+(p + q)) 
+ ^ / dp'T''(l)-^{p')T'^(l){p' - p - g)j I - 4 Tr / dpdqdfp2{q + r){p^\q + r) 



-P2(g + r))| >.(g + r)^A^(g)Ai;(r)0+(p)0(p - q - r)[T^ 



+ (g + r),A^:(p)A^(g)A,^(r) [{(0+(p + g + r)T" + T^-p - g - r)) 
+ ^ / #>+(p')T"0(p' - p - g - r)l [T^ r 



(C.2) 
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+TXfc-p))(0+(A;)-0(-A;)) + ^ / dqe<f)+{q)T^<f){q + k - p) 



x{<j)+{k)-<j){-k)) + ^ J dq{q'-{k-qY){<P^{p-k)T'^ 
+TXA; - k) + ^ J dp'dq{q^ - {k - qf)(t>^{p')T" 



X(/.(p' + A;-p)0+(g)0(g-A;)J + j rfpt/#^A;,A^(p)A^(g) - A;)r'^ 

+T"0(fc-p))(</.+ (fc + g)T^ + T^0(-fc-g)) + ^ j dp' {(j)-^ (p - k)T'' 
+T''^{k - p))(l)+{p')T^(t){p' -k-q) + ^ J dp'dcf -p + q') 

x0+(p')T^0(p' - A; - g)] + j dp{p^ - {k + pf) + 
x(0+(A;)-0(-fc)) + ^ j dq{q^-{q-kf)cP+{p)<P{k + p)<P^{q)<P{q-k) 
+k^A^^{k) [±A;2T"(0+(A;)-0(-A;)) + ^ j dp{{k + pf - p')T^<P+{k + p) 
x0(p)j + / #A;^A:,A^(A;)A^(p)[|T"(0+(A;+p)T^ + rX-A;-p)) 

(C.3) 



+^ / dqT"(j)+{q)T''(j){q~k~p) 



The divergent cutoff function p^^ enters in botli tfie expressions fIC.ip , flC.2p and flC.3l) . Tfie 
sliift of tlie momentum variable of tlie field is in general not allowed in some terms of the 
expressions (IC.ll) . flC.2l) or flC.3l) . if these terms become not regularized. For example, let us 
consider the expression fIC.ip . 

If we perform the shift of variables p p + k and g ^ g — in the term proportional 
to ^A^(p)Aj^(g), its UV behavior becomes worse. This is because, the divergent function 
p^'^{k), is now only multiplied by the product of the fields A'^^p + k)A'^{q — k), and not by 
any smeared field function of the variable k which adds the multiplicative factor p\{k). 
Having in mind the rules outlined before, one can convince oneself that the expressions ( IC.ll) . 
flC.2l) and fIC.SP are indeed regularized if the propagators of the gluons and of the ip fields 
fl2.32p are rapidly decreasing functions in momentum space. Notice that the contributions 
to the action ( 12.27^ of the terms in the second brace of ( 1C.2P can be regularized without any 
constraints on the gluons and ip propagators, knowing that the product of cutoff functions 
P2{Pi^ — P2) is a rapidly decreasing function of k'^ in momentum Euclidean space. 

The contribution to the action (12.271) of its pure YM part fl2.28p and of the first, third and 
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last terms of (12.301) . are respectively given in momentum space by the following expressions, 



S- 



YM 



I dk pf{k){k'g^,-k^kM^^{k)A''^{-k) 
+ J dkdp p^\k)p2{k)k''A>^^{k)A^p{p)A,^{-k - p) 

-j^C^p^,C^s^> I dkdpdq pl{k)A^{p)A''^{k - p)A^^{q)A^s{-k - q) (C.4) 



GG 



^ Tr y dkdpdq pl{k)i^ - {pqg,, - - k)T'^ 

+T''ct){k-p)){ct)'^{k + q)T^ + T^cl>{-k~q)) + ^ j dp\<P+ {p ~ k)T^ 
+T''^{k - p))(P+{p')TP(P{p' -k~q) + ^ [ dp'dq'(l)+{p')T''(t){p' + k - p) 



x0+(g')T^0(g' -k-q)\+j dfp,A^^{p)Af,{q)A^,{r) - k)T 

+T''<P{k - p))(0+(A; + q + r) [T^ T^] + [T^, T^]<P{-k - q - r)) 



+ |i / f/p'(0+(p-A;)T" + TXA;-p))0+(p')[T^T^]0(p'-A; 



q — r) 



+ f j dp'<P+{p')T^cP{p' -p-k){<p-^{q-k + r)[TP,T^] + [TP,T^] 
X(P{k - q - r)) + § j dp'dq'(f)+{p')T''(P{k + p' - p)<p+{q') 
x[T^T^]0(g'-A;-g-r)] + j dfd7sA''^{p)A-f,{q)A,^{r)A^s{s) 

+ [T^ T^]ct){-k - r - s)) + ^ j dp'{(f)+{p + q - k)T''T^ 

+T''T^(j){k -p- q))(j)+{p')[T^, T^W -k-r-s) + ^ j dp' dq' {p') 

xT''T^<i){k + P-P- q)<i)'^{q')[T\ T^]<i)'^{q' - k - r - s)]) (C.5) 
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^Tr J dkdppi\k)p2{k)!^-{kpg,,-k,pMi^{p)A''^{k) [i(0+(fc+p)T" 



dqk^A^,(p)A'^p(q) 



+T'-ct>{-k-p))TP + ^ j dqct>-^{k+p + q)T''ct>{q)Tf' 
A,,{k) [i(0+(p + q + A;)[r", T^] + [T", T^]ct>{-k -p- q))T^ 
+^2 J dp'(P+{p') [T«, T^](j){p' ~k-p- g)TT] I + ^ Tr y dpdqdfpl{q + r) 

p,A^:(p)A^(?)A^^(r) + g + r)T" + - g - r)) 

+^ / dp'cP+{p')T"cf>{p' -p-q-r)] [T^, T^] + / ds>l^(p)>l^(g) A^^(r) 



1 j.+ t 



ip + q + r + s)T'^T^ + T'^T^(i){-p - q - r - s)) 



j dpf<l>+(p')T-T'<l>(p' -p-q-r -s)] [T^ T^] 



(C.6) 



^IVy dkdpdq P2{k){p-,\k) - - qpkM^{p - k)T- 

+T"(t){k - p))(f>-^{q)(f>{q -k) + ^ J dp\kpp'^ - pp'k^)4)+{q)T'^4){k + q - p) 

X(/.+ (p')0(p' - k)] + {kpg,. - k.pMMA'f^il) - k)T^ + T'^Hk - p)) 

x{(P+{k + q)T^ + T'^(P{-k-q)) + ^ J dp' {p - k)T'^ + (P{k - p)) 

X(j)+{p')T^(j){p' -k-q) + ^ J dp'(j)-^{p')T''(i){k + p' + q)T'^ 

+Tl'cf>{-k - g)) + ^ I dp'dq'cf>^ip')T-cf>ik +p'- p)cf>^iq')TPcf>iq' - k - q)] 

I dp'{k,pl - k.p'^) + q - k)T<^TP + T'^T^cf^ik - p - q)) 

X(/.+ (p')0(p' -k) + ^ j dq'(j>^{q')T''T'^cP{k + q' - p - g)0+(p')0(p' - k) 

(0+(A; + r)T^ + T^<f){-k - r)) + ^ J dp'{cl>+{p + q - k)[T'', T^] + [T", T^] 
(f>{k-p-q))(f>+{p')T^(f>{p' -r-k) + ^ J dp'(f>+{p')[T'',Tf^](f>{k-p + p' -q) 
x{(f)+{k + r)T^ + T^<j){-k-r)) + ^ J dp'dq'cf)+{p')[T",T'^](j){k - p + p' - q) 



X 



X 



X(f)+{q')T'^(f){q' -k-r) 



(C.7) 
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Since the cutoff function p2 and the product of cutoff functions p2Pi^ are all rapidly decreas- 
ing functions of k"^ in momentum Euclidean space, the contributions to the action (12.271) of 
the terms (1C.5I) . (1C.6I) and (1C.7I) are readily seen to be all regularized. 

In summary, all the contributions to the action described in this appendix are regularized, 
if the cutoff function entering the smeared field f l2.32p behaves at least as fl2.33p . and if both 
the propagators of the gluons and the (p fields are rapidly decreasing functions in Euclidean 
momentum space in the IR and UV domains. Moreover the shift of momenta of the smeared 
field are allowed in all terms contributing to the action (12.271) . apart in some terms of the 
expressions ( ICUll . ((Cj) and (ICl3|l . 



Appendix D 

After integration over the boson fields y?, only the proper diagrams with internal bosons 
lines must be considered. Regarding the structure of the action S^eg (14. ip in momentum 
space, the vertices entering in such proper diagrams F can be classified in two sets. 

Let T^X^ be the set of vertices given by the parts flCUl) . fICT . flCT . flCl5l) . fIClBl) and 
flC.7p of the action (14. ip . In these vertices the derivative act in general on the boson fields 
(f. If A^^* is the number of gluon lines attached to the vertex F^*^, the order of the derivative 
acting on ip, is by inspection at most 4 — A^^* . Notice that there are no derivatives acting on 
the (p fields in a set of vertices given by the parts ( IC.SP and ( IC.SP of the action (14. ip . Thus if 
N^^ is the number of (p fields attached to the vertex F^*^^, by power counting the contribution 
of such a vertex to the UV behavior of the proper diagram F is 

r^V-A^-^i-lA-i)^^', (D.l) 

where A is the UV cutoff scale. The part S"^ (I3.12p of the action (14.10 gives the second set 
F^* of vertices. If N"^ is the number of p} fields attached to the vertex F^*, the contribution 
of this vertex to the UV behavior of F is 

Fj ~ M\K-Y'-'-\ (D.2) 

Now suppose that the proper vertex F is composed of L loops, with m vertices of type 1, n 
vertices of type 2 and I^p internal lines. Since in the UV domain the propagators of the fields 
p behave like (I3.16p . in this domain the proper vertex F scales like 

where Si^ and S'2„ are respectively a subset of m and n indices referring to the vertices 
of type 1 and 2. The total number V of vertices and the total number I^p of internal lines 
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pertaining to the proper vertex F (]D.3P are respectively given by 



V = m + n, (D.4) 

and 

Knowing that the number L of loops is related to V and [TS] by 

L = I^-V + 1, (D.6) 



the expression fID.Sp can be reduced to 

^~(^)'^"'^A-^ (D.7) 

Here the integer k is given in terms of the total number = YlieSu external gluon 

lines pertaining to the proper vertex F flD.3l) by 

k = NA-^, (D.8) 



and the dimension of F flD.7p is exactly that of A The values of k can be discussed 



according to th e val ues of Na- 

If Na > 4, (1D.8I1 shows that k is strictly positive. Then, by the dimensional argument, 



when all the external momenta scale with A, the proper vertex F (ID .70 is suppressed by 
a power of A. If Na < 4, only the two, three and four effective gluons vertices must be 
considered. At a fixed loo p ord er, t hese e ffective vertices are a sum of proper vertices given 
by ( ID. 71) . Therefore, from ( ID. 71) and ( ID.SI) . the polarization operator, whose dimension is that 

of a mass squared, seems to behave at worst like "A^ in the UV domain. However 

from the dimensional ground, Lorentz invariance and gauge invariance of the theory, this 
tensorial operator is transverse and then behaves in the UV domain like {jp) '^(9(A°). This 
means that for the two gluons effective vertex, the minimal value of k is not —2, but zero. The 
behavior of the three gluons effective vertex, is at first sight (;^) ^A. Since the dimension of 
this vectorial operator is that of a mass, the Lorentz invariance imposes that this dimension 
is that of the external gluon momenta. In that case the minimal value of k is in fa ct zero. 
For the four gluons effective vertex under consideration, the minimal value of k ( ID.SI) is zero. 
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